Double-precision calculation is used to make sure the precision of the outcome. The second order aweather format is used to the scatter of the equation. SIMPLE arithmetic which is classical in CFD calculation is used to coupling of pressure-velocity. 在此基础上采用牛顿-拉格朗日插值法进行计算,用双精度计算保证计算结果的精度,方程中各项的离散采用二阶迎风格式,压力-速度的耦合采用CFD中经典的SIMPLE算法求解。